The Viscosity of Quark-Gluon Plasma at RHIC and the LHC 
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Abstract. The specific shear viscosity (ti/>s)qgp of quark-gluon plasma (QGP) can be extracted from elliptic flow data in 
heavy-ion collisions by comparing them with the dynamical model VISHNU which couples a viscous fluid dynamic description 
of the QGP with a microscopic kinetic description of the late hadronic rescattering and freeze-out stage. A robust method for 
fixing (t] /^)qgp from the collision centrality dependence of the eccentricity-scaled charged hadron elliptic flow is presented. 
The systematic uncertainties associated with this extraction method are discussed, with specific attention to our presently 
restricted knowledge of initial conditions. With the (T}/i')QGP extracted in this way, VISHNU yields an excellent description 
of all soft-hadron data from Au-l-Au collisions at top RHIC energy. Extrapolations to Pb-l-Pb collisions at the LHC, using both 
a purely hydrodynamic approach and VISHNU, are presented and compared with recent experimental results from the ALICE 
Collaboration. The LHC data are again well described by VISHNU, with the same {rj /i)QGP value as at RHIC energies. 
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HOW TO MEASURE (t] /s)qgp 

Relativistic heavy-ion collisions create spatially deformed fireballs of hot, dense matter - in both non-central and 
(due to event-by-event shape fluctuations) central collisions. Hydrodynamics converts this initial spatial deformation 

into final state momentum anisotropics. Viscosity degrades the conversion efficiency er = |^;rj~^ ~^ ^p — fpTpfyvj 
of the fluid {x and y are the directions transverse to the beam direction z); for given initial fireball ellipticity Ex, 
the viscous suppression of the dynamically generated total momentum anisotropy Ep is monotonically related to the 
specific shear viscosity Tj /s. The observable most directly related to Ep is the total charged hadron elliptic flow [1]. 
Its distribution in pr depends on the chemical composition and pr-spectm of the various hadron species; the latter 
evolve in the hadronic stage due to continuously increasing radial flow (and so does V2{pt)), even if (as expected at 
top LHC energy [2]) Ep fully saturates in the QGP phase. When (as it happens at RHIC energies) Ep does not reach 
saturation before hadronization, dissipative hadronic dynamics [3] affects not only the distribution of Ep over hadron 
species and pr, but even the final value of Ep itself, and thus vf' from which we want to extract Tj /s. To isolate the 
QGP viscosity (77 /5)qgp we therefore need a hybrid code that couples viscous hydrodynamics of the QGP to a realistic 
model of the late hadronic stage, such as UrQMD [4], that describes its dynamics microscopically. VISHNU [5], a hybrid 
of VISH2-I-1 (Viscous Israel-Stewart Hydrodynamics in 2+1 dimensions [6]) and UrQMD, is such a code. 



(nA)QGpATRHIC 

The left panel in Fig. 1 shows that such an approach yields a universal dependence of the ellipticity-scaled total charged 
hadron elliptic flow, v*/ Ex, on the charged hadron multiplicity density per overlap area, ( 1 /S) {dNch/dy), that depends 
only on {rj /^)qgp but not on the details of the initial state model that provides Ex and S [7]. Pre-equilibrium flow and 
bulk viscous effects on these curves are small [7]. 

The QGP viscosity can be extracted from experimental data by comparing them with these universal curves. 
The right panels of Fig. 1 show this for MC-Glauber and MC-KLN initial state models (please see [7] and references 
therein for a description of these models). In both cases the slope of the data [8] is correctly reproduced; this is not the 
case for ideal nor for viscous hydrodynamics with constant T] /s. Due to the ^20% larger ellipticity of the MC-KLN 
fireballs, the magnitude of v*g^p/ex differs between the two models. Consequently, the value of ('7/5)qgp extracted 
from this comparison changes by more than a factor 2 between them. Relative to the initial fireball ellipticity all other 
model uncertainties are negligible. Without constraining Ex more precisely, (77 /s)qqp cannot be determined to better 
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FIGURE 1. (Color online) Centrality dependence of eccentricity-scaled elliptic flow [7]. 



than a factor 2 from elliptic flow data alone, irrespective of any other model improvements. Taking the MC-Glauber 
and MC-KLN models to represent a reasonable range of initial ellipticities. Fig. 1 gives 1 <47r(T7/.s')QGp <2.5 for 
temperatures Tc < T < 27^^ probed at RHIC. 

All calculations in Fig. 1 and following below were done in "single-shot" mode, where the ensemble of fluctuating 
Monte Carlo initial states was first averaged in the participant plane [7] to obtain a smooth average initial density profile 
and then evolved just once through the hydrodynamic stage. Event-by-event evolution of each fluctuating initial state 
separately and performing the ensemble average only at the end may produce somewhat less elliptic flow and thus 
slightly reduce the (i7/>s')qgp values extracted from comparison with the data [13]. The magnitude of this reduction 
depends on (tj /5)qgp [10] but is not expected to exceed (0.2-0. 3)/4;r [13]. 
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FIGURE 2. (Color online) Transverse momentum spectra (left panel) and differential elliptic flow V2(pt) (right panel) for 
identified pions and protons from 200 A GeV Au-l-Au collisions at RHIC for different centralities. Experimental data are from STAR 
and PHENIX, theoretical lines from VISHNU (see [14] for details and references). After appropriate adjustment of initial conditions 
[14] the p7-spectra are seen to be insensitive to the QGP viscosity whereas the elliptic flow depends strongly on it. For MC-Glauber 
initial conditions (upper right panels) (Tj/i)QGp = 0.08 works well, for MC-KLN (bottom right panels) (Tj/i)QGp = 0.16 works well 
for all collision centralities. In each case, changing (tj /^)qgp by 0.08 destroys the agreement between theory and data. 

VISHNU with (tj /s)qgp = for MC-Glauber and ^ for MC-KLN provides an excellent description of all aspects 
of soft (p7-<1.5GeV) hadron production (p^-spectra and differential V2(pt) for all charged hadrons together as 
well as for individual identified species) in 200 A GeV AuH-Au collisions at all but the most peripheral collision 
centrahties [14]. As an example we show in Fig. 2 pr-spectra and differential elliptic flow for identified pions and 



' It has been suggested [9, 10, 11, 12] that the ambiguity between the MC-Glauber and MC-KLN ellipticities which Ues at the origin of this 
uncertainty can be resolved by simultaneously analyzing elliptic and triangular flow, i'2 and V3 . 



protons (resonance decay contributions are included). Such a level of theoretical control is unprecedented.^ 



(tj/s)qgp ATTHE LHC 

The successful comprehensive fit of soft hadron spectra and elliptic flow in Au+Au collisions at RHIC shown in Fig. 2 
and elaborated on in more detail in Refs. [14, 15] allows for tightly constrained LHC predictions. Fig. 3 shows such 
predictions for both pure viscous hydrodynamics VISH2+1 [15] and VISHNU [16]. A straightforward extrapolation with 
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FIGURE 3. (Color online) Total charged hadron elliptic flow as function of centrality (VISHNU, left [16]) and differential elliptic 
flow for identified hadrons for 20-30% centrality (VISH2+1, right [15]) for 200AGeV Au+Au collisions at RHIC and 2.76ATeV 
Pb+Pb collisions at the LHC. Experimental data are from [17]. 

fixed (J7/i)QGP overpredicts the LHC values by 10-15%; a slight increase of {t]/s)qgp from 0.16 to 0.20 (for MC- 
KLN) gives better agreement with the ALICE data [17]. However, at LHC energies v'2 becomes sensitive to details of 
the initial shear stress profile [15], and no firm conclusion can be drawn yet whether the QGP turns more viscous (i.e. 
less strongly coupled) at higher temperatures. Furthermore, ALICE [18] has noted a discrepancy between the p/n^ 
ratio measured in Pb-nPb collisions at the LHC and the value observed by STAR in AuH-Au collisions at RHIC. The 
latter has a strong influence on the value of the chemical decoupling temperature implemented in the model. We use 
7chem = 165MeV which nicely fits the normalization of the proton spectra from STAR but overpredicts those from 
PHENIX by a factor ^ 1.5—2 (left panel in Fig. 2). The p/n^ ratio measured by ALICE at the LHC agrees with the 
PHENIX value measured at RHIC (see Fig. 7 in [18]) but is smaller by a factor ^ 1.5—2 than what is implemented 
in the LHC predictions from Refs. [15, 16]. Correspondingly, our predictions of the p spectra for PbH-Pb@LHC [15] 
overpredict the measured spectra by this factor [18]. Reducing the p/n^ ratio to the measured value will reduce the 
charged hadron elliptic flow [19]. To to the larger radial flow, this reductions is stronger at the LHC than at RHIC. This 
may account for the ^ 15% overprediction of Vj'' for (rj /i)QGp =0. 16 at the LHC seen in the left panel of Fig. 3. 

The right panel of this figure shows that, at fixed < 1 GeV, V2{pt) increases from RHIC to LHC for pions but 
decreases for all heavier hadrons. The similarity at RHIC and LHC of Vj^{pt) for the sum of all charged hadrons 
noted in Ref. [17] thus appears accidental. As a result of this shift of the elliptic flow to larger pr for heavier particles, 
which is caused by the stronger radial flow at the LHC, the mass-splitting between the viipr) curves for different mass 
hadrons grows from RHIC to LHC. This predicted growth has been confirmed by ALICE (see Fig. 6 in [20]). 

As mentioned in footnote 2, the purely hydrodynamic simulations based on VISH2+1 with constant 77/5 = 0.2 fail 
to correctly reproduce the centrality dependence of the proton elliptic flow v^ipT)- Especially in central collisions, 

(pr) is overpredicted at small pj (see Fig. 2 in [21]), i.e. the radial flow pushing the elliptic flow to higher pj (and 



^ We note that the purely hydrodynamic model VISH2+1 does almost equally well, with (T)/i)QGP = 0.2 for MC-KLN initial conditions [15], except 
for the centrality dependence of the differential elliptic flow V2(pt) for protons. We will see a similar failure of VISH2+1 for Pb+Pb collisions at 
the LHC further below. The main difference to VISHNU is that, in order to generate enough radial flow at freeze-out, VISH2+1 must be started 
earlier (To = 0.6 instead of 1.05 fm/c) because it lacks the highly dissipative hadronic phase that generates additional radial flow in VISHNU at late 
times (in the VISH2+1 simulations r}/s is held constant at 0.2 until hadronic freeze-out). The variation with collision centrality of the final balance 
between radial and elliptic flow turns out to be correct in VISHNU (where more of the radial flow develops later) but incorrect in VISH2+1 (where 
more of it is created early). 



thereby reducing V2{pt) at low pr) as generated by the model is not strong enough in central collisions. The same 
problem is seen in Fig. 4 for the extrapolation of VISH2+1 to LHC energies: whereas the pion and kaon elliptic flows 
at pj- < 1 .5 GeV/c are well described at all collision centralities, there appears to be a lack of radial flow in central and 
semi-central collisions such that the proton V2 is not pushed towards larger pr as strongly as seen in the data. This 
problem disappears in the more peripheral bins, indicating an incorrect centraUty dependence of the balance between 
radial and elliptic flow in the VISH2+1 model. 
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FIGURE 4. (Color online) Differential elliptic flow, V2{2}{pt), for pions, kaons and anti-protons from 2.76 A TeV Pb+Pb 
collisions at different centralities, as measured by ALICE (preliminary data reported at Quark Matter 2011 [20, 21]), compared 
with VISH2+1 calculations using MC-KLN initial conditions with rj /s = 0.2 [15]. 

As was the case at RHIC energies, this problem is removed when replacing the hydrodynamic description of the late 
hadronic phase with low 77/5 = 0.2 by a microscopic kinetic description using VISHNU (see Fig. 5). The calculations 
shown in Fig. 5 were done primarily to understand systematic differences between the predictions for spectra and 
elliptic flow from VISH2-I-1 and VISHNU. For this reason they were performed with the same value rj/s — 0.2 in the 
QGP phase, even though the VISHNU calculations for identified hadrons from AuH-Au at RHIC [14] had shown a slight 
preference for the smaller value {rj /.s')qgp =0. 16. Fig. 5 shows that the problem with the lack of radial flow in central 
collisions that was seen in Fig. 4 has been resolved: for the 5%-lQ% and 10%-20% centrality bins, VISHNU describes 
the differential V2 up to pj —IGeW/c almost perfectly, for all three particle species. Looking closely, one observes 
a slight (6%) underprediction of V2(pr) far all three particle species. This underprediction gets stronger in more 
peripheral collisions (reaching 9% at 50%-60% centrality) - a clear sign that we have slightly overestimated (tj /j')qgp- 
(Since, for fixed Tj /s, viscous effects increase in inverse proportion to the fireball size [6], an overestimate of Tj/i will 
lead to an underprediction of vjipr) that grows with the impact parameter of the collisions.) We are confident that, 
after reducing (tj /.?)qgp to 0. 16, the data will be well described at all collision centralities. Corresponding simulations 
are in progress. 



CONCLUSIONS 

The hybrid model VISHNU, which describes the evolution of the dense and strongly coupled quark-gluon plasma 
phase macroscopically using viscous fluid dynamics and that of the dilute late hadronic rescattering and freeze-out 
stage microscopically using a kinetic approach, provides a comprehensive quantitative description of the bulk matter 
created in relativistic heavy-ion collisions at RHIC and LHC. Transverse momentum spectra and elliptic flow of 
soft charged hadrons, pions, kaons, and protons are well reproduced at all collision centralities, with a QGP shear 



viscosity (77 /5)qgp = ^ = 0.16 if MC-KLN initial conditions are used. So far the data yield no evidence for a change 
of (t7 /s)qgp between RHIC and LHC that would reflect the different temperature ranges probed. Overall, the QGP 
liquid created in heavy-ion collisions at the LHC appears to be as strongly coupled as at RHIC energies. 




FIGURE 5. (Color online) Same preliminary data from ALICE [20, 21] as in Fig. 4, but now compared with VISHNU 
calculations with (t] /5)qgp =0.2, using the same MC-KLN initial conditions as in Fig. 3. Shown is the eccentricity-scaled elliptic 
flow, i.e. V2{2}/ev{2} for the experimental data and {vi) / {£x) for the theoretical curves. 
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